Efficiency and Fluctuation in Tight-Coupfing Model of Molecular Motor 
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A simple tight-coupling model of a molecular chemical engine is proposed. The efficiency of the 
chemical engine and its average velocity can be explicitly calculated. The diffusion constant is eval- 
uated approximately using the fluctuation theorem. Langevin simulations with stochastic boundary 
conditions are performed and the numerical results are compared with theoretical calculations. 

PACS numbers: 



Molecular motors such as myosin, kinesin, and Fi-ATPase are highly efficient for converting chemical energy into 
mechanical work, although the mechanism of the energy transduction is not fully understood. P, 0, 0, 3 Since the 
molecular motors are very small, they work under strong thermal fluctuations. Various ratchet models have been 
intensively studied, in which the fluctuations play an important role for the directional motion. [^0,0 In particular, 
Fi-ATPase is a reversible chemical engine, in which the torque and the rotary motion are induced by the chemical 
reaction: ATP^ ADP-I-P, and ATP is produced from ADP-I-P by the reverse torque. It is known that the efficiency of 
the energy conversion is close to 1. Some ratchet-type models were proposed for the rotary chemical engine. ^8. 9. lof 
We proposed a simple Langevin model with stochastic boundary conditions for the Feynman ratchet [ij and showed 
that entropy production satisfies the fluctuation theorem. In this letter, we will generalize the Langevin model 
with stochastic boundary conditions for a chemical engine such as the Fi-ATPase. The model is simple such that the 
average rotational velocity and the efficiency of energy conversion can be obtained explicitly; further, the diffusion 
constant for the Brownian motion is approximately evaluated. 

We propose a simple model of a rotary molecular motor such as Fi-ATPase under thermal fluctuation at temperature 
T . The angle of the rotary machine is denoted as x. The average velocity of the Brownian motion of the molecular 
machine is assumed to be vi in the region x < Xi and x > X2 and V2 in the region xi < x < X2. Here, we have assumed 
two special angles xi and X2. The chemical energy of ATP serves as the energy source of the molecular motor. ATP 
is adsorbed only at the angle xi. If chemical energy is provided to the motor, the motor can rotate further. A similar 
type of local reaction was also assumed in ref. 10. The adsorption of ATP is assumed to occur in a stochastic manner. 
We introduce a probability pi, which is proportional to the adsorption probability of ATP ai x — xi. The Brownian 
particle can pass through xi in the -|-a;-direction with the passing probability pi, although the Brownian particle can 
pass through xi freely in the — x-direction. This type of stochastic boundary condition was used in a model for the 
Feynman ratchet. |ll| If ATP adsorption does not occur, the Brownian particle stays within the region x < xi. The 
adsorbed ATP changes into ADP-I-P in the region xi < x < X2, and ADP is desorbed at another specific angle X2- 
The Brownian particle can pass through X2 freely in the +x direction, and the reverse motion occurs only with a 
probability p2 , which is proportional to the adsorption probability of ADP at x = X2. The average velocity after the 
release of ADP is again vi . 

The Langevin equation for x is given by 

dx 

— = ui -f £,{t), for < x < xi, a;2 < x < xo, 
dx 

= V2+^{t), fOTXi<X<X2, (1) 

where £,{t) denotes the Gaussian white noise satisfying {^{t)^{t')) = 2TS{t — t'), and the spatial period of xq = 27r/3 
is assumed, because Fi-ATPase has a rotational symmetry of the angle 2tt/S. Equation (1) can be rewritten using a 
potential as 

where U{x) — —vix for < a; < xi, U{x) ~ ~V2{x — X2) + Uq for xi < x < X2, and U{x) = — wi(a; — xo) for X2 < x < xq. 
The parameter values of xi = 7r/9 and X2 = Stt/Q are used in this paper. The probabilities pi and p2 are assumed 
to be pi = exp(— Ai?/T) exp(— A/i/T) cxp{(— wi/i — V2l2)lT} and p2 — exp{— A£^/T), where li = xq -\- Xi — X2 and 
I2 — X2 ^ xi. Figure 1(a) displays the potential U{x) for vi — V2 — 25 and AE = 20, and the constant Uq is assumed 
to be Uq = AE + Vi{xo — X2)- The unit of energy is assumed to be pN-nm, and the unit of force is pN, although 
our parameter values are not entirely realistic. The potential gap AU{x2) at x = X2 is then AE, and the potential 
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FIG. 1; (a) Potential U{x) for «i = «2 = 25 and AE = 20. (b) Time evolutions of x{t) by Eq. (3) at F = 18 and F = 24 for 
= «2 = 25, AE = 20, Afi = 50, and T = 4 



gap AU{xi) ai X = xi is AE + vih + V2l2- The probabilities pi and p2 are related to the potential gaps as follows 
Pi = exp(— A[/(a;i)/T) exp(— A/i/T) and p2 ~ exp(— AC/(x2)/T). The energy gap AU{x2) = AE can be considered 
as the activation energy for the adsorption of ADP. The energy gap AU{xi) is considered as the activation energy 
for the adsorption of ATP at x = xi. The quantity A/x denotes a difference in the chemical potential of ATP and 
ADP+P. If ATP and ADP+P are in a state of chemical equilibrium, the difference in their chemical potentials A^ is 
zero, and the total system comprising ATP,ADP, P, and the molecular machine is in a state of thermal equilibrium. 
If the concentration of ATP is greater than that in chemical equilibrium, A/i becomes negative, and the adsorption 
probability pi increases; the molecular motor then tends to rotate in the +a::-direction. The difference in the chemical 
potential is expected to obey 

A/i = Afio + Tln[ADP][P]/[ATP] 

for an ideal solution. In this case, the adsorption probability pi cx exp(— A/z/T) at x = xi is proportional to the 
concentration [ATP] of ATP. 

When the molecular motor works against the external torque F, the model equation (1) is replaced with 

dx 

— — Vi — F + ^(t), for < X < Xi, X2 < X < Xq, 
at 

dx 

— = V2- F + ^(i), for xi < X < X2. (3) 
at 

In this loaded system, if a; is increased from x = to x = xq, the potential increases from C/(0) to U{xo) = U{0) + Fxo. 
When Fxq becomes equal to the difference in the chemical potential — A/z, the chemical energy of ATP is completely 
used for the work of the uphill motion. This corresponds to the thermal equilibrium condition. The equilibrium 
condition does not depend on AE, vi, and V2 

We have performed numerical simulations of the Langevin equation (3) with the stochastic boundary conditions. 
The Euler method with a time step of 5 x 10^® was used for the numerical simulations. Figure 1(b) displays the two 
time evolutions of x{t)/xo at F = 18 and 24 for vi = V2 — 25, T — A, AE = 20, and A/i = —50. Hereafter, the 
parameter values of T = 4, AE = 20 and A/i = —50 are fixed. The Brownian motion exhibits a stepwise motion 
because the Brownian particle remains within the region nxo < x < nxo + x\ or nxo + 2:2 < x < (n + l)xo (n is an 
integer) for a long duration and occasionally jumps from uxq + x\ to nxo + X2 in a relatively short time. The average 
velocity is definitely positive for F — 18, however, it is nearly for F = 24 because F = 24 is close to the equilibrium 
condition F ^ F^ ^ -A^i/xq = 23.873. 

The corresponding Fokker-Planck equation for the probability density function P{x, t) is given as 

dP _ d{{-dU{x)/dx-F)P} d^P 

dt dx dx^' ^ ' 

The stationary probability distribution satisfying the stochastic boundary conditions: P{xi — 0) ■ pi = P{;xi + 0) and 
P{x2 - 0) = P(a;2 + 0) ■ p2 is given by 



P{x) = P(a;2 + 0) + 



c \ c 
— exp{v'i(x + Xq ~ X2)/T} -, for < a; < xi. 
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FIG. 2: Stationary probability distribution function P{x) at _F = 18 and T = 4 for (a) vi — V2 = 25 and (b) vi — 15 and 
V2 = 30. The numerical results are shown by the solid curves and the theoretical lines are shown by the dashed curves. 



— , for xi < a; < X2, 

Q 

— , for 2:2 < a; < lo, (5) 
^1 

where v'^ ^ vi — F, V2 ^ V2 — F, and the three unknown constants P{xi + 0), P{x2 + 0), and c are determined by 
the stochastic boundary conditions and the normalization condition J^" P{x)dx = 1. The constant c is proportional 
to the probability flow, and c becomes in the equilibrium state. The average velocity (v) of the Brownian motion is 
given by — c • xq. The explicit form of c is rather complicated, but it is expHcitly given by c = l/(ci + C2) where 

^ {1 - exp(t;^/i/r)}/^; + exp(A^/T){exp(t;^^i/r) - expjv'^h/T + v[h/T)}/v!, 
~ exp{(-A/i - Fxa)/T} - 1 



{p{xi + 0) + exp{«^(a; - xi)/T} 
(p{x2 + 0) + ^) exp{w;(a; - X2)/T} 



X r[{l - exp(-w'i?i/r)}/w; + exp(-A^/r) exp{(-AS - v^h - t'2/2)/T}{exp(v^?2/r) - 1}^, }] 

C2 = T{1 - exp(-w;Zi/T)}/wf - h/v[ + T{exp(v^;2/r) - l}/v'i - h/v'^. 

Figure 2(a) and 2(b) display the stationary distribution of P{x) for (a) vi — V2 — 25 and (b) vi = 15 and V2 = 30 
at _F = 18. The solid lines represent the histograms obtained from the Langevin simulations and the dashed curves 
represent eq. (5). The probability in the region xi < a; < 2:2 is very small. This is because the passing probabilities 
Pi and p2 are rather small, which is related to the stepwise motion shown in Fig. 1(b). The probability distribution 
function P{x) is an increasing function of x in the region x < xi and a:; > 0:2 for vi = V2 ~ 25 because vi > F. On 
the other hand, P{x) is a decreasing function of x for vi — 15 and V2 = 30 because vi < F. However, the average 
velocity of the Brownian motion is positive in both the cases. 

Figure 3(a) displays the average velocity (v) in two cases of vi — V2 — 25, and vi — 15, V2 = 30 as a function of F. 
The symbols represent numerical results of the Langevin simulations, and the curves represent the theoretical values. 
The average velocity becomes at _F = Fc = 50/a::o in both the cases, however, the average velocities take different 
values for F ^ Fc- The average velocity is a nonlinear function of F . The efhciency r of the energy conversion is 
defined as the ratio of the chemical potential and the work per cycle: r = Fxo/(— A/i) for {v) > 0. If the average 
velocity is negative, the external torque induces a reverse rotation. In the case of the reverse rotation, when x 
decreases from x — X2 to x — xi, ADP+P changes to ATP. This is because our molecular motor is a tight-couphng 
model, and the mechanical work is used for the production of ATP. The efficiency of the ATP production is defined 
as r = {—Afi)/Fxo for (v) < 0. Figure 3(b) represents r as a function of F. The efficiency attains the maximum 
value of r = 1 at F = Fc. The symbols denote numerical results for the two cases vi — V2 — 25 and t^i = 15, ti2 = 30. 
In the numerical simulations, r is calculated as the long time average of the ratio of Fx{t) to — A/i x AV, where Nr 
is the total rotation number of the Brownian motion. 

The fluctuation of the molecular motor as shown in Fig. 1(b) can be characterized by the diffusion constant D of 
the Brownian motion. It is dcflned as 

D = {{x(t + to) - x{to) - {x{t + to) - x{to)))^)/{2t), 

where the average is taken as the ensemble average. The diffusion constant of the Brownian motion in tilted periodic 
potentials was studied in refs. 13 and 14. Here, we evaluate the diffusion constant using the fluctuation theorem. [l5l | 
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FIG. 3: (a) Average velocity (v) of the rotary Brownian motion as a function of F for vi = V2 — 25 (denoted by rhombi) and 
for vi — 15 and V2 = 30 (denoted by crosses). The theoretical values are shown by the solid and dashed curves, (b) Efficiency 
r of the rotary molecular motor as a function of F for vi — V2 = 25 and for vi — 15 (denoted by rhombi) and V2 = 30 (denoted 
by crosses). 




FIG. 4: Diffusion constant D of the Brownian motion as a function of F for (a) vi — V2 — 25 and (b) vi = 15, V2 = 30. The 
dashed curves are the lines generated by using eq. (7). 



The fluctuation theorem is satisfied by a rather large class of stochastic systems, including our model. The probability 
P{cr) of the entropy production (t(t) during a large time interval r satisfies the relation 

Because the Brownian motion exhibits a stepwise motion, we consider a spatially@discretized lattice system with 
lattice constant xq. The Brownian motion is approximated by a random walk on the lattice. We can consider that the 
Brownian particle stays at the nth lattice site with probability P„ defined as P„ = l!ixt^^^° t)dx, where P(x, t) is 
the continuous probability distribution function. The time is also discretized as ijv — NAt with the time step At, which 
is assumed to be sufficiently small. The transition probability from the nth site to the n ± 1th site in the short time At 
is defined as p± . The probability that the Brownian particle stays at the nth site is given by po = ^ — P+ — P-- This 
random walk problem can be explicitly solved. The probability that the number of forward steps is iV+ , the number of 
backward steps is N- , and the number of times that the random walker stays at the same site is A'o = — — in 
the total steps of the random walk is given by P{N^) = Nl/ {N^IN_INq\)p^^ p^^ p^" . The entropy production in 
this process is expressed as ct = {-A^ - Fxo){N+ - N_)/T . Therefore, P(cr) Nl/ {N+\N_\No\)p+*p'^- Po° , and the 
probability of the reverse process is expressed as P{—a) = Nl/ {N+lN^lNol)p^^ p^^ Pq° . The value lnP(cr)/P(— cr) is 
calculated as (A^+ — iV_) lnp+/p_. Only if = e^^'^^"^^"^/^, the fluctuation theorem is exactly satisfied even 

for a small step number N. Therefore, the ratio p+/p- needs to satisfy p+/p- = e^~^'^~^^°'>^'^ . 

The average velocity and the diffusion constant of this random walk is therefore evaluated as 

(v) = xo{p+ -p-)/At, 

D = xl{{p++p-)-{p+-p^f}/{2At). (6) 
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For a sufRciently small Af, and p_ are also sufficiently small; therefore, the term — in the expression of 
D is negligible. In the limit of continuous time At 0, Z) is expressed as 

r, ^ ^liP+ +P-) ^ ^o('")(P+ +P-) ^ ^ot^ 

2At 2(p+-p_) 2tanh{(-A/x-Fxo)/(2r)}' 

Near F = Fc = —A^/xo, v is written as v = ^{F^ — F), where 7 = —dv/dF at F = Fg. The diffusion constant D is 
then expressed as = 7 ■ T. This is the Einstein relation at the equilibrium condition. The Einstein relation in which 
the diffusion constant is expressed by the response coefficient 7 is satisfied only near F = F^. Equation (7) is interpreted 
as a generalized relation between the diffusion constant D and the generalized response function 7(F) = v/{Fc — F). 
Figure 4 shows the numerically obtained values of D and the relation of eq. (7) as a function of F for (a) vi = V2 = 25 
and (b) Vi = 15, V2 = 30. Equation (7) is a fairly good approximation of the diffusion constant. This equation implies 
that an effective temperature can be expressed as Te// — Dl^{F) = z^/tanh(j//T), where v = xo{Fc — F)/2. This is 
a simple relation, although it is based on an approximation by the random walk. The effective temperature does not 
depend on detailed parameters such as Vi, V2 and so on. The effective temperature T^ff is greater than T, and the 
equahty is satisfied only at F = Fc, because Tg// is expressed as Te// = {(j^/r)/tanh(i^/T)} • T and 2;/tanh2; > 1 for 
any z. 

To summarize, we have proposed a simple tight-coupling model for a rotary molecular motor. Because it is a 
piecewise linear and one-dimensional model, various quantities can be calculated explicitly. Its efficiency is determined 

only by the ratio of the difference in the chemical potential — A/i and the work F -xq. Wc have calculated the diffusion 
constant approximately by using the fluctuation theorem; and the numerical results obtained are fairly well consistent 
with the approximation. 
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